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ABSTRACT 


Studies  on  cyclic  plasticity  during  the  second  year  of 
the  contract  has  developed  along  a  theoretical  approach  to 
the  problems  involved.  It  is  hoped,  however,  that  this 
approach  will  furnish  a  good  basis  for  useful  and  meaningful 
tests  on  components  operating  under  cyclic  plasticity  conditions. 

^  The  investigations  of  the  modes  of  cyclic  plastic  behaviour 
in  tubes  under  internal  pressure  and  cyclic  thermal  gradient 
have  been  continued. -^Tjhe  results  obtained  so  far  are  summarized 
in  Part  lv*>An  approximate  uniaxial  model  of  the  biaxial  tube 
problem  is  used  for  studying  the  influence  of  the  hardening  rule 
on  the  cyclic  plastic  behaviour.  Comparisons  are  made  between 
kinematic  and  isotropic  hardening  rules.  These  two  theories 
represent  the  extremes  of  cyclic  plastic  behaviour.  To  avoid 
the  drawbacks  of  the  uniaxial  model,  a  two-dimensional  numerical 
solution  utilizing  kinematic  hardening  has  been  also  developed. 
Results  obtained  are  in  better  correlation  with  experiments. 
Present  efforts  on  the  tube  problem  will  make  use  of  a  plasticity 
model  which  can  faithfully  represent  the  observed  cyclic  harden¬ 
ing  and  cyclic  creep. 

In  Part  2, 'the  cyclic  plastic  behaviour  of  specimens  used 
in  the  so-called' pulley  test  are  studied  using  linear  kinematic 
hardening  theory.  The  pulley  test  simulates  the  cyclic  thermal 
stresses  by  means  of  cyclic  bending  of  a  thin  specimen  around 
the  circumference  of  a  freely  rotating  pulley.  Theoretical 
results  for  stress  and  strain  behaviour  are  obtained  for  both 
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been  used  with  great  success  in  the  analysis  of  creeping  struc¬ 
tures  -  is  applied  to  the  analysis  of  components  operating  under 
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between  the  behaviour  of  the  structure  at  a  given  point  to  the 


result  of  a  uniaxial  test  conducted  at  the  reference  condition 
of  stress  or  strain  amplitudes.  Accordingly,  the  analysis  is 
reduced  to  the  determination  of  the  reference  test  conditions. 
This  elegant  procedure  obviates  the  need  for  extensive  testing 
programmes  aimed  at  describing  accurately  the  behaviour  of 
metals.  The  technique  is  illustrated  for  beams  under  both 
uniform  and  non-uniform  bending  and  comparisons  with  published 
test  data  show  excellent  agreement. 
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PART  1 


ANALYTICAL  AND  NUMERICAL  INVESTIGATION 
OF  THIN  TUBES  UNDER  INTERNAL  PRESSURE 
AND  CYCLIC  RADIAL  TEMPERATURE  GRADIENT 


Abstract 

The  elasto-plastic  behaviour  of  thin  tubes  under  internal 
pressure  and  cyclic  radial  thermal  gradient  is  investigated  for 
non-linear  temperature  distribution  across  the  tube  wall.  The 
following  simplifying  assumptions  are  introduced:  1)  an  equivalent 
uniaxial  model,  which  ignores  the  effect  of  axial  stress  on  the 
tube  behaviour,  is  adopted,  2)  the  r.on-linear  temperature 
distribution  is  approximated  by  means  of  a  bilinear  relation, 

3)  the  plastic  behaviour  of  the  material  is  represented  by  a 
linear  kinematic  hardening  model.  The  numerical  procedure  is 
modified  to  handle  the  isotropic  hardening  rule  in  a  similar 
manner  to  that  for  handling  kinematic  hardening. 

Results  of  this  investigation  show  that  the  bilinear 
temperature  distribution  introduces  new  modes  of  cyclic  plastic 
deformation  which  are  not  present  for  linear  distribution. 

Also,  the  bilinear  distribution  is  found  to  cause  more  stringent 
ratchetting  limits,  larger  ratchet  strains  and  larger  amplitudes 
of  cyclic  plasticity  when  compared  with  those  due  to  the  linear 
distribution. 

The  numerical  progranme  is  further  improved  so  as  to  solve 
the  tube  problem  in  a  generalized  plane  strain  form,  i.e.  in  two 
dimensions,  and  thus  the  need  for  an  equivalent  uniaxial  model 
no  longer  exists.  This  two  dimensional  numerical  solution  is 
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also  modified  to  include  isotropic  hardening.  Outputs  of 
the  programme  for  the  uniaxial  model  under  both  kinematic 
and  isotropic  hardening  rules  as  well  as  those  for  the 
biaxial  model  under  the  kinematic  hardening  rule  are  in  good 
agreement  with  analytical  and  experimental  results.  The 
biaxial  model  under  the  isotropic  hardening  rule  is  still 
under  convergence  and  accuracy  checks. 

The  results  are  also  compared  with  experimental  results 
available  in  the  literature  for  the  plastic  ratchetting  of 
thin  tubes  tested  in  a  sodium  loop  and  operating  above  the 
shakedown  limit. 

1.  Introduction 

A  unique  high  temperature  structural  design  problem  in 
the  Liquid -Metal-Fast-Breeder-Reactor  (LMFBR)  results  from 
the  frequent  thermal  transients  that  can  occur  in  heat 
exchanger  tubes  during  the  reactor  startups  and  shutdowns. 

These  intermittent  thermal  shocks  can  produce  progressive 
inelastic  deformation  (ratchetting)  and  significant  fatigue 
damage  . 

Ratchetting  of  pressurized  thin  tubes  under  cyclic  thermal 
transients  has  been  investigated  by  Miller  /l/  and  Bree  /2/ 
using  perfect  plasticity  assumptions.  Mulcahy  /3/  realized 
the  unrealistic  representation  of  the  behaviour  of  austenitic 
steels  by  perfect  plasticity  and  therefore  analysed  the  tube 
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problem  using  a  linear  kinematic  hardening  model.  These  studies 
consider  a  linear  temperature  distribution  across  the  tube  wall. 
In  practice,  however,  temperature  is  distributed  non-linearly 
across  the  wall  as  seen  in  Fig.  la  which  reproduces  the  basic 
features  of  Corum  et  al  /4/  tests  on  ratchetting  of  304  stain¬ 
less  steel  tubes.  Although  the  kinematic  hardening  model 
accounts  for  an  idealized  Bauschinger  effect,  it  does  not 
reproduce  some  equally  important  facts  of  the  cyclic  plastic 
behaviour  of  metals  such  as  cyclic  hardening,  cyclic  creep 
(material  ratchetting)  and  cyclic  relaxation. 

In  the  present  work,  the  non-linear  temperature  distribu¬ 
tion  across  the  tube  wall  is  taken  into  account  by  means  of  a 
bilinear  approximation.  A  linear  kinematic  hardening  model  is 
used  for  the  analysis  and  both  the  kinematic  and  isotropic 
hardening  models  are  employed  in  the  numerical  solution.  In 
coitmon  with  previous  work,  /l/,  /2/  and  /3/,  the  two-dimensional 
tube  problem  is  reduced  to  an  equivalent  uniaxial  one  by  ignoring 
the  effects  of  the  axial  membrane  stress  and  magnifying  the 
thermo -elastic  stress. 

The  elasto-plastic  solution  is  obtained  in  a  semi -analy¬ 
tical  form  and  the  computer  is  employed  to  perform  some  of  the 
tedious  manipulations  associated  with  the  ratchetting  regimes. 

The  results  show  that  the  bilinear  temperature  distribution 
introduces  new  modes  of  deformation  not  associated  with  the 
linear  distribution  /1,2,3/.  For  example,  operating  conditions 
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above  shakedown  will  exhibit  cyclic  plasticity  which  is 
confined  to  the  inner  skin  of  the  tube  only  as  compared 
with  cyclic  plasticity  at  both  the  inner  and  outer  skins 
for  the  linear  temperature  distribution.  Also,  the  bilinear 
distribution  introduces  more  stringent  ratchetting  limits, 
larger  ratchet  strains  and  larger  amplitudes  of  cyclic 
plasticity  at  the  inner  skin  of  the  tube  wall,  which  in  turn 
enhances  fatigue  damage  in  the  tube  /5/.  The  previous 
arguments  are  found  to  be  true  for  isotropic  hardening  as 
could  be  seen  from  the  numerical  results. 

The  two-dimensional  kinematic  hardening  solution 
technique  gave  similar  regimes  of  deformation  to  those  obtained 
from  the  uni-axial  model.  Although  the  stress  distribution 
at  a  certain  cycle  is  not  very  much  different  in  magnitudes, 
its  form  is  considerably  different  from  those  obtained  for 
the  uniaxial  model.  This  may  be  attributed  to  the  bi-axiality 
effects.  Also  the  ratchet  strains,  cumulative  strains  and 
amplitudes  of  cyclic  plasticity  are  markedly  lower  than  those 
for  the  equivalent  uniaxial  model. 

The  present  elasto-plastic  results  are  used  to  analyse 
the  tube  ratchetting  tests  carried  out  by  Corum  et  al  /4/  on 
thin  tubes  made  from  initially  annealed  304  stainless  steel 
under  severe  thermal  downshocks  between  1100°F  and  800°F, 

Fig.  la.  In  these  tests,  each  thermal  shock  is  followed  by  a  hold 
period  at  1100°F  which  causes  accumulation  of  time-dependent 
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creep  strains.  This  creep  strains  can  be  evaluated  using 
the  technique  developed  by  Ainsworth  /6.7/  which  yields  an 
upper  bound  on  creep  deformation  per  cycle  when  load  variar:o 
exceed  the  shakedown  limit . 


2 .  The  Elasto-Plastic  Problem 

Consider  a  thin  tube  (inner  radius  a  and  outer  radius  b 
under  an  internal  pressure  P  ami  a  cyclic  thermal  gradient  of 
maximum  value  AT  distributed  non-lineuriy  across  its  wall 
thickness  t.  The  non-linear  temperatui e  distribution  may  be 
approximated  by  means  of  a  bilinear  relation  such  that  : 

T  (x)  =  0  forOsx^h 

T{x)  =  -  AT  (x-h)/(l-h)  for  h  $  x  *  1 

where  x  =  (b-r)/t,  i.e.  x  =  0  at  r  =  b  and  x  -  l  at  >:  -  a 
Elastic  thermal  stresses  are  given  by  (/«/'): 


_  _  .  Ea  A  T 

°e  -  °z  “  ~  (1-h)2TI^v~ 


for  0  i  x  h 


(2) 


„  .  ,  En  A  T  ,  X-h  EaAT.  .  .  „  .  , 

°e  az  ~  (3_h)2(l-v)  1-h  1-v  h  1 


These  stresses  are  cyclically  superimposed  upon  the  membrane 
stresses  due  to  internal  pressure,  which  are  given  by 


_  _  _  PD 

°e  -  2oz  "  2t 


(4) 


where  D  is  the  mean  diameter  of  the  tube. 

As  argued  by  Bree  /2/,  the  effect  of  oz  is  opposite 
that  of  o0,  and  hence  ignoring  a  should  yield  conservative 
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strain  estimates.  This  reduces  the  two  dimensional  problem 
to  that  of  a  slab  prevented  from  bending  and  subjected  to  a 
steady  mechanical  stress  =  PD/2t  in  addition  to  the 
temperature  cycle  given  by  equations  (1),  Fig.  lb.  The 
corresponding  uniaxial  thermal  stresses  are, 

at(x)  =  -  (l-h)at  for  o  <  x  <  h 

(5) 

a  (x)  =  -  (l-h)o  +  2  (x-h)  h  <  x  <  1 

C  11  ( 1 — h )  °t 

where  EaAT/(l-v)  is  one  half  of  the  total  thermal 

stress  across  the  tube  wall.  Fig.  lb.  An  operating  condition 

will  therefore  be  defined  in  terms  of  ap  and  0fc. 

The  elasto-plastic  analysis  is  carried  out  using  a 
linear  kinematic  hardening  material  model,  Fiq.  2,  where  the 
slope  of  the  plastic  portion  is  assumed  to  be  a  constant  BE. 
This  yields  the  flow  rule  as, 

%  -  w-  6°  -  I d-  <6> 

where  K  ■  (l-B)/8  is  the  material  hardening  parameter.  For 
simplicity,  all  material  constants  E,  v,  a,  K  and  are 
assumed  to  be  temperature  independent. 

In  the  isothermal  condition,  the  uniaxial  model  is 
subjected  to  a  uniform  stress  o^,  and  upon  application  of  the 
thermal  gradient  during  the  firt  half  cycle,  the  stresses  will 
increase  in  the  cold  fibers  and  decrease  in  the  hot  ones. 
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Consequently,  two  distinct  stress  fields  are  possible,  the 
first  involves  tensile  yielding  only.  Fig.  3a,  and  the  second 
involves  both  tensile  and  compressive  yielding,  Fib.  3b.  The 
elasto-plastic  stress-strain  relations  for  the  stress  field 
of  Fig.  3a  are: 


el  =  ai<x)/E  for  °<x<h 

G1  =  "  2ofc(x-h)/(l-h)]/E  for  h<x<y1  (7) 

e 1  =  [a1(x)  -  2at(x-h)/(l-h)  +  Kfo^x)  -  a  )]  for  y1<x<l 

where  y^  is  the  plastic  front.  Equilibrium  and  compatibility 
conditions  are  expressed  by: 


o  = 


o 


/  aL(x)  dx 


and 


e  ^  =  constant 


(8) 


Since  a^ty.^)  =  ay,  the  constant  c1  becomes  equal  to  [ay-2o+ (yl-h)  / 
(l-h)]/E  and  the  integrals  in  the  equilibrium  equation  can  now  be 
evaluated  to  yield  y^  as: 

yl  *  "  F  +  F  ^(1+K) (1+Kh2)  +  K(l+K) (1-h) (l-ap/oy) .a  /afc  (9) 

Fully  elastic  behaviour  is  assured  if  y^  z  1  which  implies  (using 
Eq.  (9) )  that: 

o t  *  (o  -  o  )/(l+h)  (10) 

Compressive  yielding  occurs  at  x  =  o  if  0^(0)  <  -  which  is 
equivalent  to  (using  Eq.  (9)): 

(°t  /  °y)  2  +  [(1+K)  d-op/oy)  “  2  ( 1+Kh)  ]  .  (ofc  /  0y)  /  (1-h)  -  K  >  0 


(ID 
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Operating  conditions  which  satisfy  inequality  (11)  will  cause 
both  tensile  and  compressive  yielding  during  the  first  half 
cycle  (Fig.  3b)  and  the  solution  of  the  corresponding  elasto- 
plastic  problem  yields  the  following  expression  for  the  strain 

el: 

G1  a|>y- 2ot(y1-h)/(l-h)]/E  =  -[ay+2ot(z1-h)/(l-h)]/E,  where 
(1+h2)  (a./  a  )  2+  (1-h)  (1+K)  (l-o  /a  )  .  a  Jo  +K  (1-h2) 

y  =  - l—X. - E— * - 5— X -  (12) 

2.ot/o  (0t/a  +  K(l-hJ) 

Subsequent  thermal  cycling  gives  rise  to  either  elastic  shake- 
down,  cyclic  plasticity  or  ratchetting,  as  shown  in  Fig.  (4). 

3 .  Elastic  Shakedown,  Cyclic  Plasticity  and  Ratchetting 

For  operating  conditions  which  satisfy  Eq.  (11) ,  the 
stress  fields  shown  in  Figs.  5a  and  5b  are  possible  during  the 
second  half  cycle.  The  stress  field  of  Fig.  5a  involves  cyclic 
plasticity  at  the  inner  skin  of  the  tube  while  the  stress  field 
of  Fig.  5b  involves  cyclic  plasticity  at  both  the  inner  and 
outer  skin  of  the  tube. 

The  elasto-plastic  analysis  of  the  stress  field  in  Fig.  5a 
yields  the  plastic  front  y 2  given  by  Eq.  (13) 

y2  =  "  k  +  i  /(1+K) (1+Kh?)+2K(1+K) (1-h) 0/0 t  (13) 

Shakedown  will  occur  if  y2  *.  1  which  is  satisfied  if 

ot  <  2  0  /(1+h)  (14) 

and  such  shakedown  behaviour  is  termed  as  (Fig.  4) .  In 

Fig.  5a,  cyclic  plasticity  will  not  occur  at  the  outer  skin 

(x=o)  ifAa-(o)  <  2a  and  such  condition  is  satisfied  if 
*  Y 


L  , 
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<Tt  <  <ry  [/ (K-2Kh-l)  2  +  4K(l-h)2  -  (K  -  2Kh-l )  ]  /  ( 1-h ,  (15' 

Thus,  thermal  stresses  which  satisfy  Eq.  (15)  but  viulii. 
F.q .  (14)  produce  cyclic  plasticity  at  the  ii.ik;  skin  oniy,  and 
such  behaviour  is  called  in  Fig.  4.  Note  that  for  h  =  o, 
Eqs.  (14)  and  (15)  become  identical,  viz.  ^  2 and  regime 
is  not  present  as  in  Mulcahy's  analysis  /3/  tor  linear 
temperature  distribution. 

The  stress  field  of  Fig.  5b  will  occur  when  thermal 
stresses  violate  Eq.  (15)  and  cyclic  plasticity  takes  place 
at  both  the  inner  and  outer  skin  of  the  tube,  Such  behaviour 
is  termed  F2  in  Fig.  4.  It  can  be  shown  that  F2  will  continue 
during  subsequent  temperature  cycling  if 

[<l+Kh2)-(l4K)ap/oy|  .  (ot/oy)2tK(l-h)  [3-2 (1  +  K; op/oy] .afc/ay+2K2 

(1-h)  £  o  (16) 

as  shown  in  Fig.  4. 

Cyclic  plasticity  may  also  take  place  when  the  response 
of  the  model  exhibits  tensile  yielding  only  during  the  first 
half  cycle  (Fig.  3a).  A  possible  second  half  cycle  stress 
field  which  corresponds  to  such  behaviour  is  shown  in  Fig.  5c 
and  the  corresponding  plastic  front  y2  is  given  by  Eq.  17. 

y2  “  “  i  +  K  ^<1+K> Yl+Kh5)  +  2K(K+1)  (1-h)  .oy/ot 


k 


(17) 
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This  mode  of  cyclic  plasticity  is  similar  to  i . 
reversed  plasticity  is  confined  to  the  inner  skin  cr.iy ,  Such 
behaviour  is  termed  in  Fig,  4  and  it  can  be  shown  that  trie 
condition  for  its  continuation  during  subsequent  thermal 
cycling  is 

[(1+K)  (l+o  /a  )  2-4  (1+Kh2) ]  .  (ofc/oy)  2-2K  (1-h)  (3-,y'oy) . a^/a  M  K“  (1-hV <1  .+K)  *c 


Shakedown  will  occur  if  v-  *  5  and  c 'o)  <  <j  ana  these 

■  x  / 


conditions  are  satisfied  if, 


at  * 


4 (1+K) 

K (1+h) (1-h^) 


(a  -a  )  and  n  s  2o  /(1  +  h) 

y  p  v 


This  shakedown  behaviour  is  termed  £?  in  Fig,  4. 


Ratchetting  takes  place  for  operating  conditions  which 
violate  Eqs.  (16),  (18)  or  (19)  as  can  be  seen  from  Fig.  4  for 
the  ratchetting  regimes  R-^ ,  R2*  R^  and  I<  .  ,  In  the  ratchettinq 
regime  R^,  which  takes  place  at  thermal  stresses  below  the 
shakedown  limit  (afc  =  2cry/(l+h)),  a  cyclic  stress  field  which 
produces  cyclic  permanent  deformation  is  shown  in  Fig.  5d . 

The  accumulated  ratchet  strain  is  obtained  thro.;-/: 
numerical  evaluation  of  the  recurrence  relations  for  the 
corresponding  yn  and  y  ^  by  a  simple  computer  algorithm. 
Elastic  shakedown  is  approached  asymptotically  as  the  cyclic 
steady  state  following  an  initial  phase  of  transient  ratchett.r. 

In  regime  transient  ratchetting  takes  place  while 

cyclic  plasticity  is  occurring  at  both  the  inner  and  outer  skin 
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of  the  tube  (Fig.  5e)  until  a  cyclic  steady  state  similar  to 
F2  is  approached  asymptotically.  In  regime  ,  the  transient 
response  consists  of  two  consecutive  phases  of  ratchetting 
mechanisms.  In  the  first  phase,  ratchetting  takes  place 
according  to  a  mechanism  similar  to  that  present  in  until  a 
sufficient  degree  of  hardening  is  attained  such  that  the 
mechanism  shown  in  Fig.  5f  is  activated.  In  this  second  phase 
of  R^,  ratchetting  is  developed  while  cyclic  plasticity  is 
taking  place  at  the  inner  skin  of  the  tube.  Clearly,  the  cyclic 

steady  state  for  will  be  similar  to  that  present  in  and  F.^. 

In  the  ratchetting  regime  R^ ,  ratchetting  is  developed  through 
three  consecutive  phases.  R^  is  employed  in  the  first  phase,  the 
second  phase  is  similar  to  the  second  phase  of  R^  and  the  third 

phase  is  the  same  as  R2.  Consequently,  the  cyclic  steady  state 

for  F2,  R2  and  R^  are  similar. 

When  a  linear  temperature  distribution  (h=o)  is  assumed 
across  the  tube  wall,  the  resulting  cyclic  steady  state  is 
either  shakedown  for  ofc  <  2oy  or  cyclic  plasticity  at-  both  the 
inner  and  outer  skin  of  the  tube  for  >  20y.  It  is  now  clear 

that  due  to  the  bilinear  temperature  distribution  (h  >  o) ,  new 
modes  of  behaviour  (F^,  F^  and  )  are  discovered  in  which  cyclic 
plasticity  is  confined  to  the  inner  skin  of  the  tube. 

4.  Discussion  of  Results 


The  results  of  previous  sections  were  assembled  into  a 
computer  program  described  in  detail  in  a  previous  report  /9/ 
whose  input  is  the  problem  parameters  1<K<®,  o<h<0.5 
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and  the  operating  condition  as  defined  by  o  <  and 

o  <  crt  <  <» .  The  program  is  provided  with  10  routines  which 

yield  stresses  and  strains  corresponding  to  E,  S^,  S2,  ,  F2, 

F3,  ,  R2>  and  R4.  In  the  ratchetting  regimes,  the  cyclic 

steady  state  is  assumed  to  be  reached  wheni£p(x=o)  does  not 
change  by  more  than  10-^  per  cycle.  Examples  of  the  transient 
ratchetting  behaviour  in  R^,  R2,  R^  and  R^  are  shown  in  Fig.  6a. 

It  is  seen  that  most  of  the  ratchet  strain  accumulates  during 
the  first  20  cycles  or  so.  The  number  of  cycles  required  to 
reach  the  steady  state  depends  mainly  on  the  hardening  parameter 
K.  As  K  becomes  larger,  the  material  plastic  behaviour  approaches 
that  of  perfect  plasticity  and  consequently  a  larger  number  of 
cycles  is  required  before  the  steady  state  is  reached.  However, 
for  austenitic  and  many  high  strength  steels  K  *  40  -  50  which 
imply  that  the  steady  state  will  be  reached  after  the  first 
20  to  30  cycles. 

In  Fig.  6b  the  effect  of  h  is  clear.  The  ratchet  strains 
as  well  as  the  cumulative  strain  are  considerably  lower  in  the 
case  of  linear  temperature  distribution  (h=o)  than  for  the 
bilinear  distribution.  This  is  true  for  both  kinematic  harden¬ 
ing  and  isotropic  hardening  results  as  could  be  seen  from  Fig.  6c. 

An  overall  view  of  the  results  may  be  gained  from  the 
contours  of  total  ratchet  strains  as  shown  in  Fig.  7  for  a 
thin  tube  made  from  304SS  whose  monotonic  properties  at  900°F 
are  used  to  calculate  pressures,  temperatures  and  strains.  It 
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may  be  noted  that  the  contours  of  equal  ratchet  strains  are 
almost  parallel  to  the  ratcheting  limit  (the  line  which 
separates  R  regions  from  S  and  F  regions)  and,  correspondingly, 
ratchet  strains  are  more  sensitive  to  changes  in  pressures  than 
changes  in  temperatures.  In  fact,  ratchet  strains  are  reduced 
slightly  as  AT  becomes  very  large  (for  AT  >  250°F  in  Fig.  7). 
However,  these  large  temperature  differences  are  not  likely 
to  be  encountered  in  practice. 

The  important  result  which  emerges  from  the  present  work 
is  the  effect  of  the  bilinear  temperature  distribution  parameter 
h  on  the  tube  behaviour.  As  stated  earlier,  for  h  >  o  new  modes 
of  deformation  (F^,  F^  and  R^)  are  introduced  in  which  cyclic 
plasticity  is  confined  to  the  inner  skin  of  the  tube.  This  is 
shown  in  Fig.  8a  where  it  can  be  seen  how  h  affects  the  form  of 
stress  distribution  for  a  certain  set  of  loading  parameters 
(those  used  in  Corum's  tests).  The  cyclic  plasticity  at  the 
inner  surface  of  the  tube  observed  for  h  =  0.328  is  not  present 
for  h  =  o.  This  is  true  for  both  kinematic  and  isotropic  harden¬ 
ing  solutions. 

Figure  8a  shows  the  effect  of  h  on  the  ratchetting  limit 
and  it  is  seen  that  ratchetting  starts  to  take  place  at  smaller 
mechanical  loads  as  h  becomes  larger.  This  effect  is  more 
prominent  above  the  shakedown  limit  (ofc  =  2a^) .  Fig.  8c  shows 
the  effect  of  h  on  asymptotic  ratchet  strain  at  =  0.4cjy, 
o ^  =  3cjy.  It  is  seen  that  ratchet  strain  attains  a  maximum 
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value  at  h  =  0.27.  Thus,  an  analysis  based  on  linear  tempera¬ 
ture  distribution  is  likely  to  underestimate  the  likelihood  of 
ratchetting  as  well  as  the  value  of  ratchet  strain.  The  effect 
of  h  on  the  amplitudes  of  cyclic  plasticity  which  are  of 
interest  when  designing  against  low  cycle  fatigue  is  shown  in 
Fig.  8c.  At  h  *  o,  cyclic  plasticity  of  equal  amplitude  takes 
place  at  both  the  inner  and  outer  skin  of  the  tube  and  as  h 
becomes  larger, cyclic  plasticity  amplitudes  increase  at  the 
inner  skin  (x=l)  and  decrease  at  the  outer  skin  (x=o) .  For 
example,  at  h  =  o,  Aep(o)  =  Ae^Cl)  =  1.5cy/E  while  at  h  =  0.3, 

Ae  (1)  =  3  o  /E,Ae  (c)  =  o  which  means  tnat  the  higher  value 

c  y  hr 

of  h  doubles  the  amplitude  of  cyclic  plasticity  at  x  =  1. 

This  can  be  shown  to  reduce  the  fatigue  life  by  about  thirty 
percent  (using  the  Cof f in-Manson  relation  /5/:Ae  .N?'5  =  constant). 

tr  ■** 

5 .  Analysis  of  Inelastic  Benchmark  Ratchetting  Tests 

Two  thermal  ratchetting  tests  have  been  conducted  by 
Corum  et  al  /A/  on  304SS  thin  tubes.  In  the  first  test,  a 
tube  was  subjected  to  13  thermal  downshocks  between  1100°F  and 
800°F  in  a  sodium  loop  under  internal  pressure  of  700  psi  which 
is  reduced  to  zero  for  2  hrs  subsequent  to  each  thermal  shock 
in  order  to  transfer  the  liquid  sodium  from  the  drain  to  the 
source  tank.  Following  each  thermal  shock,  the  tube  was  kept 
at  1100°F  for  160  hr.  The  measured  temperature  distribution 
across  the  tube  wall  (Fig.  la)  was  non-linear  and  may  be 
reasonably  approximated  by  a  bilinear  relation  with  h  =  0.328 


as  discussed  earlier.  The  maximum  temperature  difference 
across  the  wall  was  about  167°F  and  corresponds  to  an  outer 
skin  temperature  of  1000°F  and  an  inner  skin  temperature  of 
830°F. 

Using  the  monotonic  properties  at  an  average  temperature 
of  900°F,  the  full  solution  map  is  constructed  as  shown 
previously  in  Fig.  7  which  shows  that  the  operating  conditions 
of  this  test  lie  within  the  ratchetting  regime  R^.  The  cyclic 
accumulation  of  ratchet  strain  calculated  on  the  basis  of  mono¬ 
tonic  data  at  900°F  is  shown  in  Fig.  9  as  dotted  lines. 

Ainsworth  /6,7/  developed  a  method  by  which  an  upper  bound 
on  creep  strain  per  cycle  can  be  determined  for  structures 
subjected  to  load  variations  above  the  shakedown  limit.  The 
creep  strain  is  bounded  by  the  cyclic  plasticity  solution  of 
a  similar  structure  that  does  not  creep.  Thus,  considering  a 
uniaxial  tube  model  subjected  to  internal  pressure  P(t)  and 
cyclic  temperature  T(x,t),  both  of  cycle  t,  and  denoting  the 
cyclic  plasticity  solution  for  a  similar  structure,  which  does 
not  creep,  subjected  to  the  loading  conditions  P(t)  +  p  and 
T(x, t)  by  a*(x,t),  an  upper  bound  for  axial  creep  displacement 
is  given  by 

Au  <  —  /  /  D  o  (x,t) J  dV  dt,  (20} 

where  p  is  an  additional  constant  positive  load  yet  to  be 
determined  and  D  is  the  rate  of  creep  energy  dissipation. 
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The  bound  is  evaluated  at  the  value  of  p  which  minimizes  the 

right-hand  side  of  Eq.  (20)  .  Note  that  the  earlier  analysis 

★ 

of  the  elasto-plastic  problem  yields  the  stress  field  a  (x,t). 
Assuming  that  creep  is  governed  by  Norton's  law  and  substitu- 
ting  for  D  by  Ka  ,  Eq.  (20)  becomes, 


*u<  K±  (^)n+1  i  /V(x)+1dx, 


(21) 


Most  of  the  creep  strain  takes  place  at  the  hold 
temperature  of  1100°F  and  the  short  term  creep  rates  /10/  at 
1000-2000  hrs  are  used  for  evaluation  of  the  bound  for  the 
test.  The  stress  index  n  is  obtained  as  n  =  3.655,  /10/,  at 
the  level  of  stresses  present  in  the  tests. 

The  final  predictions  of  plastic  and  creep  ratchetting 
are  compared  with  experimental  results  of  the  test  in  Fig.  9. 
The  values  of  yield  stress  are  a =  11.15  ksi  during  the 
first  half  cycle  and  a ^  -  14.47  ksi  during  subsequent  cycling 
which  corresponds  to  a  cyclic  plasticity  amplitude  of  +  0.05% 
as  obtained  from  initial  calculations  based  on  the  monotonic 
data  at  900°F.  The  creep  strain  is  0.0156%  per  cycle  which 
corresponds  to  an  optimum  additional  pressure  p  =  346.4  psi. 
Fig.  9  shows  that  the  theoretical  predictions  are  in  excess 
of  the  observed  strains  by  about  40%  for  the  test. 

The  predicted  ratchet  strains  are  found  to  be  very 
large  when  compared  with  the  experimentally  observed  ratchet 
strains  (the  experimental  results  for  the  test  are  also  shown 
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in  Fig.  9).  This  overestimation  of  observed  strains  is  seen 
to  be  due  to  the  use  of  the  monotonic  stress  -  strain  data 
and  the  kinematic  hardening  rule  which  ignore  the  cyclic 
hardening  properties  of  the  tube  material. 

The  cyclic  hardening  behaviour  of  304SS  is  character¬ 
ized  in  Fig.  10  which  shows  the  percent  increase  in  yield 
stress  against  the  number  of  cycles  and  the  amplitude  of 
reversed  plastic  strain.  Yield  stress  is  defined  here  on 
the  basis  of  bilinear  representation  of  cyclic  hardening 
loops  /10/.  The  method  suggested  here  to  account  for  cyclic 
hardening  makes  use  of  the  raonotonic  yield  stress  during  the 
first  half  cycle  only  and  the  13th  cycle  yield  stress  for 
subsequent  cycling.  The  13t.h  cycle  yield  stress  is  evaluated 
at  an  amplitude  of  cyclic  plasticity  bz^  which  is  representative 
of  the  test  conditions.  An  upper  bound  for  bz^  can  be 
evaluated  on  the  basis  of  the  monotonic  data.  This  approach, 
however,  ignores  the  effects  of  other  important  phenomena 
such  as  cyclic  creep  (material  ratchetting)  and  cyclic  relaxa¬ 
tion,  which  are  present  when  plastic  cycling  takes  place  around 
a  mean  value  of  stress  and  strain  respectively. 

One  obvious  reason  behind  the  conservatism  of  the 
predictions  is  the  modelling  of  the  biaxial  tube  problem  by 
means  of  an  equivalent  uniaxial  model  in  which  the  equivalence 
is  based  on  the  stresses  given  in  Fig.  9.  The  uniaxial  model 
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is  likely  to  predict  larger  hoop  strain  accumulation.  The 
results  of  the  biaxial  kinematic  hardening  programme  are 
also  shown  in  Fig.  9.  These  seem  to  be  the  nearest  to  the 
experimental  results  though  creep  strains  are  not  taken  into 
consideration.  This  may  be  attributed  to  the  lack  of  inter¬ 
action  between  plasticity  and  creep  in  the  classical  theories 
of  inelastic  deformations  used  in  the  present  analysis. 

It  is  also  possible  that  creep  strains  are  reduced  due 
to  prior  cyclic  hardening  /ll/.  However,  in  the  absence  of 
reliable  constitutive  relations  which  accounts  for  this  inter¬ 
action,  classical  theories  provide  conservative  strain 
estimates  for  the  conditions  of  the  present  tests. 

6.  Conclusions 

The  main  conclusions  of  the  present  work  can  be  summarized 
as  follows  : 

1.  Comparisons  between  the  modes  of  behaviour  for  linear  and 
bilinear  temperature  distributions  across  the  tube  wall 
show  that  a  bilinear  distribution  introduces  new  modes  of 
deformation  in  which  cyclic  plasticity  is  confined  to  the 
inner  skin  of  the  tube. 

2.  Also,  bilinear  distributions  are  shown  to  introduce  more 
stringent  ratchetting  limits  above  the  shakedown  limit, 
larger  ratchet  strains  and  larger  cyclic  plasticity  ampli¬ 
tudes  at  the  inner  skin  of  the  tube  when  compared  with 
corresponding  linear  distribution. 


19 


3.  It  is  demonstrated  that  the  cyclic  hardening  phenomenon 
should  be  taken  into  account  in  the  analysis  of  ratchetting 
components  above  the  shakedown  limit.  An  approximate 
method  which  accounts  for  cyclic  hardening  is  used  in  the 
present  work. 

4.  The  equivalent  uniaxial  model  of  thin  tubes  yields  conserva¬ 
tive  estimates  of  observed  strains  in  tube  ratchetting 
tests  and  the  biaxial  model  was  found  to  be  very  close  to 
the  experimental  results,  though  cyclic  creep  and  relaxation 
are  ignored. 

5.  The  inclusion  of  the  isotropic  hardening  rule  into  the 
computer  programme  increased  the  ratchet  strains  and  changed 
the  solution  map  completely.  This  is  to  be  further  investiga¬ 
ted  in  future  work. 
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Fig. 7.  Full  solution  moo  and  contours  of  total  ratchet 

strain  for  a  tube  with  D/t«2151 . h«0  32B  using  mono- 
tonic.  proposes  of  304  SS  at  900  F  I er«11.15»1Q3  PS) . 

E s  23.3  »10^  psi.  a«11,10'6/  F  Ks38.S6,  V  »0.3) 
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PART  2 


KINEMATIC  HARDENING  ANALYSIS  OF  RATCHET 
STRAIN  IN  THE  "PULLEY  TEST" 


Abstract 

The  classical  Bree  problem-which  represents  an 
uniaxial  model  of  a  thin  tube  subjected  to  combined 
internal  pressure  and  cyclic  thermal  stress  across 
its  wall-can  be  simulated  by  means  of  the  pulley  test 
in  which  a  wire  or  strip  specimen  is  subjected  to 
combined  steady  tensile  stress  and  cyclic  bending 
stress.  In  this  paper,  accuimilation  of  ratchet  strain 
in  the  pulley  test  is  investigated  using  a  linear  kine¬ 
matic  hardening  material  model  from  which  perfect 
plasticity  can  be  generated  as  a  special  case.  'rhe 
results  of  the  investigation  show  that  asymptotic 
ratchet  strains  are  linearly  related  to  the  excess  in 
mean  stress  above  its  value  g-^*  at  the  ratchetting 
limit  regardless  of  the  thermal  stress  amplitude. 
Comparisons  with  test  results  on  copper  wire  specimens- 
which  exhibit  non-linear  hardening  rate-confirm  the 
qualitative  validity  of  this  simple  relation.  Divia- 
tions  between  theory  and  experiment  are  attributed  to 
metallic  cyclic  creep.  Further,  perfect  plasticity 
results  are  shown  to  be  well  predicted  by  a  linearized 
lower  bound  estimate. 

1.  INTRODUCTION 

Incremental  strain  growth  (or  ratchetting)  and  low 
cycle  fatigue  damage  can  be  major  causes  of  failure  of 
mechanical  com  onents  subjected  to  cyclic  thermal  load¬ 
ing.  One  of  the  common  ratchetting  problems  in  power 
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generation  industries  is  that  of  a  thin-wa]  pire 
subjected  to  combined  internal  pressure  and  cyclic 
thermal  gradient  across  its  wall  which  is  known  as 
the  "Bree  problem"  [1].  Testing  of  such  components 
under  realistic  thermal  conditions  usually  involves 
prohibitive  cost  in  terms  of  enuinnent  and  strain 
measuring  devices  [2 ,3J  .  Such  difficulties  have  moti¬ 
vated  the  development  of  simpler  tests  in  which  cyclic  the¬ 
rmal  stresses  are  simulated  by  means  of  cyclic  bend¬ 
ing.  This  is  achieved  by  what  is  termed  as  the 
"pulley  test"  [4,5,6]  in  which  a  wire  or  a  strip  spe¬ 
cimen  of  the  material  whose  behaviour  is  to  be  inves¬ 
tigated  is  cyclically  bent  around  the  circumference  of 
a  freely  rotating  pulley  while  subjected  to  steady 
tensile  load  as  illustrated  in  Fia.  (la).  Considering 
the  segment  ABC  where  AB=BC=Tr?  and  ?  is  the  radius 
of  the  pulley,  when  the  dead  weiohts  W  are  in  position 
(1),  AB  is  subjected  to  combined  tension  and  bending 
while  BC  is  subjected  to  tension  only.  Alternatively, 
when  the  dead  weiqhts  are  moved  to  position  <  2) ,  BC 
will  be  subjected  to  combined  tension  and  bending 
while  AB  is  subjected  to  tension  only.  .In  this  manner, 
the  cyclic  vertical  movement  of  the  dead  weights  brings 
about  a  combination  of  steady  tension  and  cyclic  bending 
in  segments  AB  and  BC  of  the  specimen.  The  interaction 
between  elastic  and  plastic  strains  during  this  loading 
cycle  causes  incremental  strain  growth  which  manifiots 
itself  in  cyclic  lengthening  of  the  specimen. 

The  pulley  test  has  been  utilized  by  Meqahed  et 
al  [4]  as  a  simple  test  by  which  the  concepts  of  rat¬ 
cheting  and  shakedown  can  be  illustrated  experimentally. 
Figure  (lb)  shows  a  sample  of  the  results  obtained  for 
the  cyclic  strain  qrowth  of  copper  wire  specimens.  It 
is  seen  that  the  observed  behaviour  resembles  to  a  qreat 
extent  the  ordinary  creep  behaviour,  i.e.  an  initial 
primary  stage  followed  by  a  secondary  stage  of  constant 
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cyclic  strain  rate.  The  prediction  of  observed  steady 
state  rates  by  means  of  rapid  cycle  solutions  has  been 
discussed  by  Meqahed  et  al  15],  Ng  161  conducted 
ratchetting  tests  on  thin  strip  sppecimens  using  a 
pulley  apparatus  to  investigate  the  behaviour  of  316 
stainless  steel  at  room  and  elevated  temperatures. 

In  this  paper,  theoretical  solutions  for  the  deve¬ 
lopment  of  primary  ratchet  strain  in  both  wire  and 
strip  specimens  are  generated  using  the  kinematic 
hardening  plasticity  theory.  This  theory  has  the  ad¬ 
vantage  of  simulating  the  Bauschinger  effect  observed 
in  many  metals  but  lacks  the  capability  of  reproducing 
the  cyclic  phenomena  of  metals,  viz,  cyclic  hardening, 
cyclic  creep  and  cyclic  relaxation  which  have  a  pro¬ 
nounced  influence  on  the  behaviour  of  s  ructural  com¬ 
ponents  subjected  to  cyclic  loading  17].  However,  the 
solutions  presented  here  provide  an  estimate  of  cyclic 
strain  growth  as  predicted  by  a  classical  theory  of 
plasticity  commonly  employed  in  many  inelastic  finite 
element  packages  [8]. 

2 .  THE  EQUIVALENT  BREE  PROBLEM 

Elastic  bending  stresses  which  arise  in  a  specimen 
of  thickness  d  bent  around  a  pulley  of  radius  ?  are 
given  by  +  <rt  where  cr^  =  Ed/2?  .  The  stress  varia¬ 
tion  along  the  thickness  of  the  specimen  is  given  bv: 


tr(x)  =  crp  +  o-t  ( 2x-l ) 


(1) 


where  ,7.^  is  the  direct  tensile  stress  due  to  dead 
weights  and  x  is  a  dimensionless  position  variable 
along  the  thickness  d  with  x  =  0  at  the  point  of  speci¬ 
men  contact  with  the  pulley  and  x  =  1  at  the  most  outer 
fiber  (Fig.  2a).  It  is  clear  therefore  that  the  cyclic 
vertical  travel  of  the  dead  weights  gives  rise  to  a 
cyclic  elastic  stress  field  which  varies  between 


36 


cr(x)  =  cr  and  c r(x)  =  cr  +  cr  (2x-l)  .  It  can  be 
P  p  t 

shown  that  an  equivalent  problem-which  produces  the 
same  cyclic  stress  field  as  above-is  that  of  an 
assembly  of  infinite  number  of  bars  which  are  allowed 
to  deform  in  the  axial  direction  only  under  the  com¬ 
bined  effect  of  a  uniform  tensile  stress  cr  and  a 

th  ^ 

cyclic  thermal  strain  £  (x)  which  cyclically  varies 

between  -  2i .  x/E  and  zero  (Fig. 2b)  .  In  the  multi¬ 
bar  assembly,  the  constant  total  strain  is  the  sum  of 
elastic  and  thermal  strain  components,  i.e.  £(x)  =  tr(x) 
/E-2  crt.  x/E  =  constant  =  C  whivh  gives  the  stress 
history  cr(x)  =  EC+2c^.x.  The  equilibrium  conditions 

is  „jtr(x)  dA  =  cr  .  A  where  A  is  the  cross-sectional 
A-'  p  o  o 

area  of  the  specimens  and  dA  is  the  area  of  a  bar  at 

the  relative  position  x.  Note  that  A  =  1,  dA  =  dx 

for  a  strip  specimen  of  unit  width  and  Aq  =  nV4,dA=2 

] x  ( 1-x)  .  dx  for  a  wire  specimen  of  unit  diameter 

(Fig. 2c).  The  constant  total  strain  C  can  be  obtained 

by  satisfying  the  above  condition  of  equilibrium  and 

the  elastic  stress  field  o r(x)  is  obtained  as  <r(x)  = 

trr  +  or.  .  (2x-l)  which  is  the  same  as  eq .  (1)  .  Upon 

P  t 

removal  of  thermal  strain,  cr(x)  becomes  equal  to  o-p  • 

3 .  MODES  OF  CYCLIC  PLASTIC  BEHAVIOUR 

In  the  elasto-plastic  solution,  a  linear  flow 
rule  and  a  kinematic  hardening  rule  are  employed.  Thus, 
plastic  strain  increments  are  given  by 

d£P  =  £  do-  (2) 

where  K  =  (1-$)/$  and  £  E  is  the  slope  of  the  linearized 
plastic  portion  of  the  stress-strain  curve.  The  kine¬ 
matic  hardening  rule  implies  that  the  size  of  the 
yield  surface  remains  constant  at  2 oy  while  its  center 
changes  continuously  due  to  plastic  deformation. 
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Upon  the  application  of  thermal  strain  £th(x) 
for  the  first  time  (first  cooling),  the  resulting 
elastoplastic  stress  field  in  the  multi-bar  assembly 
may  involve  tensile  yielding  at  the  outer  fibers 
(Fig. 3a)  or  both  tensile  yileding  at  the  outer  fibers 
and  compressive  yielding  at  the  inner  fibers  (Fig. 3b). 
For  the  case  shown  in  Fig.  (3a) ,  the  stress-strain 
relations  are  given  by: 

Ej_  =  (cr L  ( x )  -2trt.x]/E  for  0<x  <  ax 

(3) 

l1  =  (tri(x)  -2o-fc.x+k  (o^  (x) -o-  )  )  /E  for  a^x^l 

where  a^  is  the  elasto-plastic  interface.  3ince 
°l^al^  =  ^y'  t'_ie  total  strain  is  therefore  given 
by : 


Ex  =  [o^  -  2crt.  a^  J  /E  (4) 

substit  .ting  in  eq .  (3)  provides  the  stress  field 

^(x)  which  upon  insertion  in  the  equilibrium  condi¬ 
tion  and  evaluating  the  resulting  integrals  yields 
the  following  expressions  for  a^: 


16K 

IT 


X<V 


(K+l) 


^y  ~  ^p  _  K+2 
°t  " 


(2ai-l) 


(5a) 


for  the  wire  specimens  where  the  function  %(a)  is 
defined  by: 


% (a)  =  [a(l-a) ] 3/2+16  (2a-l) [sin-1 (2a-l) +2 (2a-l) 


ya(l-a) J/3 


(6) 


and,  _ 

a^  =  [-1+  J  ( 1+k)  ( 1+k  cry/o^. .  (l-o^/irt)  ]  /K  (5b) 

for  the  strip  specimen. 
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Equations  (5a)  and  (5b)  can  be  used  to  determine 

the  operating  conditions  (tr  ,  which  give  rise  to  a 

given  state  of  deformation.  For  example  fully  elastic 

behaviour  (which  is  equivalent  to  a.^  ^  1)  will  occur 

if  ir  +  tr  <  cr  ,  which  could  be  easily  verified  from 
p  t  y 

Fig.  (2a).  Also,  compressive  yielding  will  commence 
at  x  =  0  if  tr^(O)  obtained  from  eq .  (3)  becomes  less 

than  -  it  .  This  condition  is  equivalent  to  a^cr/o^ 
which  upon  substitution  in  eqs .  (5a)  and  (5b)  yields 


<l*k)  <rp  +  <ry 


for  the  wire  specimen  and, 


2 

{-£■)  +  [  K+l )  (1- 

Y 


o-D  crt 

f2"-  <r 

y  y 


-  K^O 


for  the  strip  specimen.  Equations  (7a, b)  define  the 
locii  of  operating  conditions  corresponding  to  the 
limiting  state  between  the  two  stress  fields  shown  in 
Figs .  ( 3a ,b) . 

Analysis  of  the  stress  field  shown  in  Fig.  (3b) 
is  dealt  with  in  a  similar  manner  to  above.  The  stress- 
strain  relation  are  given  by: 

£  i  =[  (x)  -2cr  .x+Kt^  (x)+rv)  )/E  for  0  x  bx 


t1  *[«rL(x)-2«r;.x]/E 


for  b^  ^  x  <  a^ 


Ex  =[£r1(x)  -2  or  .x+K  (o^  (x)  -«^)  ] /E  for  a±  <  x  ^  1 

Noting  thatcr1(bI)  =  -  v  and  ^(a^  =cry«  follows 
that; 


Cl  -  [V2Vai)/E  = 


-  [o'y+20't.b1J/E 
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Note  that  the  width  of  the  elastic  core  a,-b,  =1  -  tr./cr,  . 

-LX  2 

The  plastic  fronts  a^  and  b.^  can  be  determined  from  the 
following  relations: 


16K 


CX<b1)  -X(ax)] 


(K+l) 


for  the  wire  specimen  and 

°rt/0y  +(K+1)(1-V+K  Y/,rt 
bl  2[V/-«r-  +  KT  - - 


for  the  strip  specimen. 


(10a) 


(10b) 


Incremental  strain  growth  or  ratcheting  will  take 
place  only  if  the  center  line  of  the  specimen  (x=0.5) 
experiences  tensile  plastic  deformation  due  to  first 
cooling.  Corresponding  operating  condition  can  there¬ 
fore  be  determined  by  imposing  the  condition  a^=0.5 
on  the  above  results.  This  procedure  defines  the 
ratchetting  limits  as: 


^  +  WT  for  't*  2<ry  (11*> 

(K+l)  _K  pl4^[X(0.5-o-y/o-t) -X(0.5) 

for  2  (lib) 


for  the  wire  specimen  and. 


°p  i  K  ^t 

o-y  T  K+T  o-y 


4  l 


°'p.0"t  1  ®*t  /  K 

- K+^  ~  - 


Y 


o-y  ~  K+ 


for  crt  4  20^  (12a) 

for  cr.  ^  2cr  (12b) 

y 


for  the  strip  specimen.  Operating  conditions  which 
satisfy  eqs .  (11a, 12a)  will  give  rise  to  elastic 
shakedown  (regions  SI,  S2)  and  those  which  satisfy 
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eqs  (11a, 12b)  give  rise  to  reversed  plasticity  '  . 

F)  as  shown  in  the  interaction  or  Bree  diagram  of 
Fig.  (4).  The  steady  tensile  stress  which  defines  the 
ratchetting  limit  as  determined  from  the  above  inequa¬ 
lities  is  denoted  by  0"  *  .  Note  that  ratchetting  limits 
corresponding  to  perfect  plasticity  material  (Fig. 4b) 
can  be  generated  from  the  above  results  by  letting 
K  —•0o.  This  yields  the  familiar  results  of  Bree  [1] 
for  the  strip  specimen  and  the  analogous  results  for 
a  circular  Bree  element.  It  is  noted  that  the  rat¬ 
chetting  limits  o-p*  corresponding  to  perfect  plasticity 
do  not  differ  considerably  from  kinematic  hardening 
solution  for  K  =  50  (Fig.  4a)  which  is  characteristic 
of  many  strain  hardening  materials. 

4 •  ACCUMULATION  OF  RATCHET  STRAIN 
a .  Kinematic  Hardening  Solution 

Elasto -plastic  analysis  of  loading  cycles  subse¬ 
quent  to  the  first  cooling  shows  that  ratchet  strain 
accumulates  by  means  of  two  distinct  mechanisms  deno¬ 
ted  as  R-,  and  R2  (Fig.  4a).  The  ratchetting  mechanism 

r  operates  at  cr  <2cr.  as  well  as  during  the  initial 
1  t  ^ 

phase  of  the  intermediate  regime  R^— *-R.,.  The  mecha¬ 
nism  P,  operates  at  as  well  as  the  second 

t  y 

phase  of  regime  R1  — *-R2.  Defining  a  typical  cycle  N 
by  the  (n-1)  th  cooling  and  nth  heating  half  cycles, 
where  n=2N.  Figure  (5a)  shows  the  cyclic  stress 
field  present  during  the  Nth  transient  cycle  and  the 
steady  state  of  the  mechanism  and  Fig.  (5b)  shows 
the  analogous  stress  field  for  mechanism  .  It  is 
seen  that  R^  will  invariably  exhibit  tensile  yielding 
which  takes  place  in  the  fibers  an_-^<x  <  1  during 
cooling  and  in  the  fibers  0<x<an  during  heating. 

The  analysis  shows  that  ^  <  0 . 5  <  during 

the  transient  stage  but  as  cycling  continues 
both  an  and  an  approach  the  value  of  0.5  at  the 
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steady  state.  Figure  (6)  summarizes  the  results  obtained  from 
the  elasto-plastic  analysis  of  the  ratchetting  mechanisms.  It 
is  shown  that  the  plastic  strain  increment  during  the  Nth 
cycle  N  is  given  by: 

SN  -  <V  »„-!>•  2VE  (13> 

Hence,  6^  -*-0  at  the  steady  state  and  the  behaviour  approaches 
an  elastic  shakedonw  state. 

In  regime  R2 ,  the  cyclic  stress  field  exhibits  both 
tensile  and  compressive  yielding  during  each  half  cycle  as 
shown  in  Fig.  (5b)  .  Again  at  the  steady  state  (a^  ***  an  j=0.5) 

,  '>0  and  the  plastic  fronts  corresponding  to  compressive 

yielding  bn  and  bn_^  approach  0.5+ir .,/ort,  0.5-oy/crt  respectively. 
Thus,  the  behaviour  at  the  stationary  state  is  that  of  a  re¬ 
versed  plasticity  at  the  inner  and  outer  fibers  of  the  assembly 
and  there  will  be  an  associated  elastic  core  of  width  2°-y/V 
The  spatial  variations  of  plastic  strain  at  the-  stationary 
state  are  illustrated  in  Fig. (5) .  It  is  clear  that  the  cyclic 
reduction  in  the  increment  of  ratchet  strain  is  due  to  strain 
hardening  and  the  presence  of  cyclic  plasticity  above  trt=  2o^ 
is  due  to  the  employment  of  the  kinematic  hardening  rule.  In 
the  intermediate  regime  R^  ->-1*2'  the  behaviour  changes  from 
R^  to  R2  when  sufficient  hardening  has  evlolved  such  that  re¬ 
verse  yielding  in  compression  starts  to  take  place. 

The  manner  by  which  the  final  results  are  presented 
in  Fig. (6)  is  very  convenient  for  computer  programming .  These 
results  have  been  used  to  determine  the  cyclic  accumulation 
of  ratchet  strain  for  operating  conditions  representative  of 
the  three  regimes;  R^ ,  R2  and  R2  as  shown  in  Fig.  (7).  For 
K  =  50 ,  which  is  characteristic  of  many  alloy  steels,  the 
behaviour  approaches  a  steady  state  in  about  50  cycles.  The 
computer  program  has  also  been  used  to  generate  the  contours 
of  the  asymptotic  values  of  ratchet  strain  G*5  (0.5)  in  the  Bree 
diagraun  as  shown  in  Fig.  (a a)  .  Further  discussions  of  ratchet 

I 
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strain  contours  and  their  relations  with  the  limit,  met-:  -,trs 
o"p*  will  be  presented  in  the  following  section  of  the  paper. 

b .  Perfect  Plasticity  Solution 

The  perfect  plasticity  solution  for  the  accumulation 
of  ratchet  strain  can  be  generated  from  the  kinematic 
hardening  results  by  letting  the  strain  hardening  co¬ 
efficient  K-*-00  .  This  has  the  following  effects? 

1.  The  intermediate  ratchet  mechanism  R^-> 
will  not  be  present  in  the  manner  described 
earlier  since  the  switch  from  to  R2  is 
primarly  due  to  strain  hardening.  Hence,  this 
intermediate  regime  will  involve  tensile  yielding 
only  in  exactly  the  same  manner  described  earlier 
for  mechanism  employed  below  er^.  —  2^  as  shown 
in  Fig.  (4b) . 

2.  Ratchet  strain  will  accumulate  in  equal  incre¬ 
ments  per  cycle  and  hence  total  ratchet  strain 

is  unbounded  qs  shown  in  the  examples  of  Fig. (7). 

Letting  K-*°°  in  the  recurrence  relations  presented  in 
the  algorithm  shown  in  Fig, (6)  yields  the  following 
closed  form  expressions  for  ratchet  strain  in  the  strip 
specimen ; 

<$N  =  2tr/E.  [l-2  )frt  ]  for  regime  (R-j^) 

2  7  (HI 

=  20-  /E.[cr  cr  /er  -1  ]  for  regime  (R0) 

N  y  p  t  y  2 

which  is  the  familiar  result  obtained  by  Bree  [1  ]. 

The  corresponding  results  for  the  wire  specimen  can  not 
be  reduced  to  closed  forms  and  will  still  have  to 
programmed . 

5.  DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

The  results  of  ratchet  strain  accumulation  obtained  in 
previous  section  are  now  compared  with  simplified  bounding 
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results  [9  ]  as  well  as  experimental  observer! .c.s  ine 

ratchetting  limits  c r  *  separate  the  ratchettLnc 

regime  from  shakedown  regime  for  cr  <  2<r  and  ratcbettinc 

regime  from  reversed  plasticity  for  2t  .  For  perfect 

plasticity  models,  Ponter  and  Cocks  f  9  j  derived  a  linearized 

relationship  between  the  increase  in  ratchet  displacement  per 

cycle  and  the  increases  in  mean  and  cyclic  thermal  stresses; 

cr  ,tr.  above  the  shakedown  limit.  This  result  is  used  to  ob- 
P 

tain  a  simple  lower  bound  which  for  kinematically  determinate 
structures  shows  that  for  moderate  thermal  loading,  the  cyclic 
increment  of  ratchet  displacement  is  greater  than  four  times 
the  elastic  displacement  of  the  body  if  it  were  subjected  only 


to  the  increase  in  mean  stress  a.bc*-e  c . 


c  cr 


+ 

:rr  ~  rr 
P  P 


This 


result  is  strictly  correct  for  loads  oelow  the  shakedown  limit 
cr^  =  2(Ty.  Ponter  and  Cocks  [  9J  compared  the  lower  bound  result 
with  exact  solution  for  the  classical  3ree  problem  which  is 
equivalent  to  the  strip  specimen  ir.  '-.he  present  work. 


Figures  (8a,b)  compare  exact  solutions  w.tn  the  lineari¬ 
zed  lower  bound  for  wire  and  strip  specimens  3-  =  cr., ,ov=2P* 

*  -  X  ^  y 

which  lie  below  the  shakedown  limit  (c r  =  2cr  )  as  well  as  at 

* 

cr  —  40-  . cr,  =  6<r  which  exceeds  cr  .  The  linearized  lower 
t  Y  t  *  y  t 

bounds  at  <rt>crt  are  calculated  by  excluding  the  regions  of 
reversed  plasticity  at  the  inner  and  outer  fibers.  For  example, 
when  a  strip  specimen  is  subjected  to  cr^  =  4c*  ,  reversed 
plasticity  takes  place  at  the  inner  fibers  0<  xx 0 . 25  and  the 
outer  fibers  0.75  <x<l  and  the  fibers  0.25\  x<0.75  experie¬ 
nce  elastic  deformation  only  (Fig. 5).  The  lower  bound  is  cal¬ 
culated  for  an  equivalent  elastic  strip  of  vidth=  0.5  and 
subjected  to  the  original  axial  force  which  means  that  the 
effective  mean  stress  is  doubled.  The  comparisons  shown  in 
Figs.  (8a,b)  indicate  that  the  lower  bounds  are  quite  close 
to  exact  sol  itions  except  whena^  approaches  cr^  . 


Investigations  of  kinematic  hardening  predictions  of 

asymptotic  ratchet  strain  at  the  core  of  the  specimen  shows 

thatE^  (0.5)  is  linearly  related  to  the  excess  in  cr  above 
*  K*h  p 

cr  in  the  following  manner; 


•*.4 


£  f  ,  O.b)  -  K/'E .  i\  cr  K/F  .  (tr  -  cr  )  ( ;  =  . 

k .  n  p  p  d 

at  any  value  of  cr_  eitner  below  or  above  tr  *  =  2<t  .  The 
dependence  of  plastic  strain  on  cr  „  is  implied  by  tne  vari¬ 
ation  of  <r  *  with  tr.  .  It  is  interesting  to  note  that  the 
p  t 

rate  of  change  of  plastic  strain  with  respect  to  cr  equals 
■  P 

the  linear  strain  hardening  coefficient  K/E  This  means 

that  the  stress  state  at  tne  core  of  the  specimens  (0 . 5 )  is 

given  by  the  stress-strain  curve  of  the  material.  In  fact, 

examination  of  the  stress  states  at  x  =  0.5  as  obtained  from 

the  computer  program  shows  that: 


°"(  0.5)  =  cr  +-<r  _  cr  *  (16') 

y  d  p 

Although  the  above  results  are  obtained  for  linear  kine¬ 
matic  strain  u-.rdcmng  material,  analogous  results  for  non¬ 
linear  hardening  can  be  derived.  This  may  be  checked  bv  con¬ 
sidering  the-  ratenetting  tests  conducted  by  Meg  ah  ad  et  al 
f4,5]  on  copter  wires  whose  operating  conditions  and  ratche- 
ttina  limits  tr*(  <=  shown  in  Fig.  (5a;.  The  stress-strain 
cui  e  of  the  copper  wires  exhibits  non-linear  hardening  which 
can  oe  fit"-  c  t  •  the  power  lav;; 


where  cr  is  the  elastic  limit,  A,  n  are  material  constants. 

p 

Hence,  the  result  forr  (0.5)  which  should  be  applied  to  non¬ 
linear  hardening  material  in  an  analogous  manner  to  eq.(15) 
should  be : 

£P(0.5)r-A  (tr  -cr*)1/n  (18) 

P  P 

recalling  the  sample  of  observed  ratchet  strain  shown  in 
Fig.(l)  and  noting  that  the  steady  ratchetting  rate  attained 
after  the  initial  primary  stage  is  mainly  due  to  cyclic  creep 
phenomena  [5]  ,  ratchet  strain^  (0.5)  appearing  in  eg. (18)  is 
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interpreted  as  the  primary  ratchet  strain  continent  ti'xy .  i  . 
Figure (9b)  shows  the  experimental  primary  ratchet  strain  at 
three  levels  of  cr^  against  ( -  <^*)  compared  with  the  stress- 
strain  curve  of  the  material.  It  is  clearly  seen  that  the 

p 

experimental  £  (0.5)  lie  on  a  single  curve  which  supports  the 


universal  dependence  of  ratchet  strain  oncr-<r*  .  However,  the 

P  P 

experimental  curve  of£p(0.5)*  lies  above  the  stress-strain 
curve  of  the  material  and  the  aiviation  is  attributed  to  pri¬ 


mary  cyclic  creep  [5}  . 
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APPLICATIONS  OF  REFERENCE  STRESS  METHOD  (RSM) 


IN  CYCLIC  PLASTICITY  ANALYSIS 


Abstract 

The  concept  of  reference  stress  method  (P.SM)  relies  on 
relating  the  inelastic  deformation  behaviour  of  a  structure 
at  a  specific  point  with  the  behaviour  of  a  uniaxial  specimen 
under  reference  test  conditions.  The  RSM  has  proved  success¬ 
ful  in  time-dependent  creep  analysis  since  it  obviates  the 
need  for  excessive  collection  of  material  data.  In  this 
paper,  the  RSM  is  applied  to  evaluate  the  plastic  behaviour 
of  components  subjected  to  cyclic  loading  conditions  where 
phenomena  such  as  cyclic  hardening,  cyclic  softening,  cyclic 
creep  and  cyclic  relaxation  are  dominant.  Several  examples 
of  beams  under  both  uniform  and  non-uniform  bending  are  dis¬ 
cussed.  Comparisons  between  the  results  of  RSM  and  several 
independent  test  results  show  very  close  agreement. 


1.  Reference  Stress  Method  in  Creep  Analysis 

The  reference  stress  method  RSM  has  been  under  develop¬ 
ment  during  the  past  twenty  years  as  an  approximate  technique 
for  analysis  of  creeping  structures.  Excellent  reviews  of  the 
RSM  may  be  found  in  the  works  of  Kraus  [1] ,  Ponter(2]and  Siml3l  . 
In  order  to  recapture  the  basic  concepts  of  the  technique,  a 
beam  with  rectangular  cross-section  (b  x  2h)  is  considered. 


The  beam  is  subjected  to  steady  bending  moment  M  and  its 
material  creeps  according  to  Norton's  creep  law? 


<r  -  A<r 


n 


(1) 


Assuming  that  plane  sections  remain  plane  during  creep,  the 
strain-curvature  relations  are; 


v  (y)  *  K.y  =  A  tr  (y) 


n 


<21 


where  K  is  the  curvature  rate  and  y  is  a  coordinate  along  the 
depth  of  the  beam.  Substituting  or  (y)  as  obtainedfrom  eq.{2) 
in  the  equilibrium  condition, 


M  *  2b  f  o-  (y)  .y.dy 
0J 


(3) 


yields  the  curvature  rate  K  as 

K  =  A  f  mT  1  fen+f" 


nn 

M 

- - - X' 

Lbh  . 


1 

’  h  * 


2n 


(4) 


Equation  (4)  can  be  rewritten  as: 


K  *  A«ron  .  i  .  f  (n) 


(3) 


where 


st  o<- 


M 


bh 


,  f(n)  * 


2ntl"Jr 
[2n-cx  J 


(6) 


The  function  f (n)  can  be  made  to  be  weakly  dependent  on  the 
creep  exponent  n  by  proper  choice  of  the  value of  c *  .  Consi¬ 
dering  that  3*5^  £  11  for  most  creeping  metals  and  making 
f  ( 3 >  *=  f(ll)  yields  o  value  of  ***  1.003.  The  corresponding 
values  of  the  function  f(n)  -  as  listed  in  Table  (1)  _ 
exhibit  very  little  dependence  on  n.  Adopting  an  average 
value  of  f  *  1.578,  the  curvature  rate  is  written  as: 


K  *  1.578  ❖[  1.003  M/bh2] 


(7) 
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Table  1  Values  of  f(n)  for  beam  under  pure  bending 


n 

f  (n) 

K  [RSM]  /K  [Exact] 

3 

1.572 

1.004 

5 

1.584 

0.996 

7 

1.584 

0.996 

9 

1.579 

0.999 

11 

1.572 

1.004 

Hence, the  beam  curvature  rate  is. obtained  directly  from  the 

results  of  a  single  creep  test  conducted  at  a  reference  stress 
2 

g-  =  1.003  M/bh  .  The  errors  envolved  in  the  approximate  solution 
are  really  quite  neglegible  as  can  be  seen  from  the  ratio 
K  [RSMj  /K  [Exact]  listed  in  Table  (1)  . 

The  identification  of  the  reference  stress  in  the  above 
procedure  relies  on  making  the  function  f (n)  weakly  dependent  on 
the  creep  exponent  n  and  hence  may  prove  difficult  for  components 
which  have  no  exact  creep  solution.  An  alternative  and  effective 
way  for  determining  o-Q  depends  on  scaling  down  the  limit  load 
solution  corresponding  to  perfectly  plastic  material  with  yield 
stress  (Ty.  It  can  be  shown  that  the  ratio  tr0/<Ty  equals  m/Ml 
where  is  the  limit  bending  moment,  i.e. 

°-0/«ry  =  M/^  (8) 

In  cases  where  no  exact  limit  load  expression  is  available , tests 
on  perfectly  plastic  models  of  the  prototype  structures  may  be 
utilized. 

The  RSM  has  been  verified  experimentally  for  beams  of  both 
rectangular  and  tubular  sections  under  both  steady  and  variable 
loads.  Creep  of  components  such  as  plates,  tubes  and  panels 
have  been  studied  using  the  RSM.  In  the  following  section, the 
RSM  is  applied  to  analysis  of  structures  made  from  cyclic  har¬ 
dening  material  under  conditions  of  repeated  load reversals. 
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2.  Cyclic  Plasticity  Analysis 

Initially  annealed  metals  are  known  to  exhibit  cyclic  har¬ 
dening  when  subjected  to  repeated  reversals  of  load  (4)  .  The 
hardening  rate  can  be  quite  large  during  the  first  few  cycles 
but  decreases  contineously  until  a  cyclic  steady  state  is 
reached  early  in  fatigue  life.  The  steady  state  of  a  cyclically 
hardening  metal  is  characterized  by  the  cyclic  stress-strain 
curve  which  connects  the  cusps  of  settelled  hysteresis  loops 
obtained  at  different  strain  (or  stress)  amplitudes.  The 
transient  cyclic  hardening  characteristics  can  be  represented 
by  a  family  of  cyclic  stress-strain  curves  which  connects  the 
cusps  of  hysteresis  loops  at  equal  number  of  cycles  and  diffe¬ 
rent  strain  amplitudes  as  shown  in  Fig.(l)  for  99.9%  pure 
copper . 

When  a  structure  made  from  cyclic  hardening  metal  is 
subjected  to  repeated  load  reversals,  its  behaviour  will  reflect 
the  hardening  characteristics  of  its  material  in  terms  of  its 
load  carrying  capacity  or  deformation  amplitudes.  In  order  to 
apply  the  reference  stress  method  to  cyclic  plasticity,  a 
simple  example  of  a  rectangular  beam  (bx2h)  subjected  to  cyclic 
curvature  *K  is  considered.  It  is  required  to  determine  the 
conditions  of  a  reference  uniaxial  test  which  can  be  used  to 
characterize  the  beam  behaviour.  The  family  of  cyclic  curves 
shown  in  Fig.  (1)  may  be  represented  by 

o-  =  B  ,  m  =  m  (N)  (9) 

where  N  is  the  number  of  cycles.  Taking  the  range  of  0.1/m 
4.0.5  which  fits  the  observed  hardening,  the  problem  can  be 
solved  in  a  manner  analogus  to  the  creep  problems.  Again 
assuming  that  plane  sections  remain  plane  and  neglecting  elastic 
strains,  the  moment  curvature  relation  at  the  instants  of  cur¬ 
vature  reversals  is  obtained  as: 

M  -  bh2  .  july  .  B  (  ±  kh)m  (10) 

which  can  be  written  as: 
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M  =  bh  .  f  (m)  .  cr(  ±  e  ) 

o 


where 


£  =«  kh  ,  f  (m)  =  - - - j,. 

°  (m  +  2  )  oc 


(11! 


Enforcement  of  f(0.1)  =  f(0.5)  yields  °<  =  0.674  and  the  corres¬ 
ponding  variation  of  f  with  m  is  as  listed  in  Table  (2) .  Hence 
the  moment  curvature  relation  can  be  written  as 


M  =  0.993  bh2  cr  (+  0.647  kh) 


(12) 


i.e.  the  bending  moment  is  determined  directly  from  the  results 

of  a  reference  uniaxial  test  conducted  at  f  =  +  0.647  times 

o  - 

the  surface  strain  (kh) . 


Table  2:  Values  of  f (m)  for  beam  under  cyclic  curvature 


m 

f  On) 

M  [RSM]  /M  [Exact] 

0.1 

0.995 

0.998 

0 . 2 

0.992 

1.0 

0.3 

0.991 

1.001 

0.4 

0.992 

0.998 

0.5 

0.995 

0.963 

Another  example  of  the  application  of  RSM  s  provided  by 
considering  the  cyclic  torsion  of  a  solid  bar  (rad.  =  a,  length=L) 
subjected  to  repeated  reversals  of  twist  angle  0.  Assuming  the 
validity  of  von-Mises  relations,  the  reference  strain  solution 
for  the  twisting  moment  T  is: 

T  =  1.205  a3  tr  (+  £  )  ,  £^  =  0.426a  0 /L  (13) 

when  the  structure  is  subjected  to  controlled  history  of  pres¬ 
cribed  deformation,  the  reference  uniaxial  test  is  conducted 
under  reversed  strain  conditions.  On  the  other  hand, when  the 
structure  is  subjected  to  controlled  history  of  prescribed 
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loads,  the  referenee  uniaxial  test  is  conducted  under  reversed 
stress  conditions.  Reference  stress  relations  in  the  latter 
case  will  be  similar  to  those  used  in  creep  analysis. 

The  above  results  describing  cyclic  pure  bending  of  beams 
are  utilized  to  generate  solutions  for  a  host  of  allied  beam 
bending  problems,  e.g.  cantilevers  and  simply  supported  beams. 
The  final  results  are  summarized  in  table  (3) .  The  range  of 
0.l4m40.5  is  used  throughout. 


Table  3.  Reference  Test  Solutions  in  Cyclic  Plasticity 


Structure  Loading  history 

Reference  stress (or  strain) solution 

Rectangular  +  K 

beam 

(bx2h)  +  M 

M=0.993bh2.  cr  [Eq]  ,  ZQ=  ±  0.647.kh 
K=1.55?i  •  £  o rQ=  +  1.005M/bh2 

Solid  bar  +  0 

(rad=a, 

Length=L)  +  T 

T=1 . 205  a3  tr[E0]  ,  EQ=  +  0.426.ae/L 

0=2 . 354L/a-  £  r  0*o=  +  0.829.T/a3 

Cantilever  +6  at  tip 

(L) under 

Tip  Load  +  P  at  tip 

2 

P=l.  158^|  .«r[E0]  ,  En=  +  1.966.£h/L2 
5=0.483—  *  c  [f"ol  '  ‘’cT  ±  0.876PL/bh2 

Cantilever  +  S  at  tip 

(L) under 

distributed  .  ,  . ,  , 

load  ±q  (unlform> 

<^=2.39^|  .0-[lo]  ,  Eq=  +  2.804.Sh/L2 

L  2 

S=0.331^-  -E(o-0]  *  ctq=  +  0.428.qL2/bh2 

Simple  Beam  +&  at  center 

(2L)  under 

central  load  +P  at  center 

2 

P=2 . 316—jr^--  r  [C0  J  ,  E0*  +  1.9665h/L2 

5  =0.483^-  •£.  jj-0]  ,  <rQ=  +  0.438.PL/bh2 

3 ■  Comparison  With  Test  Results 

The  refernece  test  results  are  compared  with  three  sets  of 
experimental  data  for  beam  bending.  The  first  set  of  tests  are 
conducted  by  Oldroyd ( 5 )  on  rectangular  copper  beams  under  pure 
bending  in  an  attempt  to  derive  the  uniaxial  stress-strain  curve 
from  the  moment-curvature  curve  by  use  of  Nadai's  bending  formula 
16  I  .  Only  monotonic  loading  is  considered  in  the  tests.  Hence, 
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a  monotonic  stress-strain  curve  and  the  central  deflection  .vs.  end 
load  in  a  four-point  beam  bending  test  are  determined  as  shown 
in  Fig.  (2).  Application  of  eq.  (13)  yields  the  reference  test 
predictions  which  are  seen  to  be  in  close  agreement  with  the 
results  of  the  bending  test. 

The  second  set  of  test  data  (7 J  involves  cyclic  bending 
of  rectangular  beams  made  from  initially  annealed  99.9%  pure 
copper  until  the  cyclic  steady  state  is  reached.  The  tests  are 
conducted  under  repeated  reversals  of  curvature  and  the  bending 
moment  is  determined.  The  reference  test  calculation  are  based 
upon  the  uniaxial  cyclic  stress-strain  curve.  An  exact  numeri¬ 
cal  solution  for  the  moment-curvature  relation  is  also  determined 
on  the  basis  of  representing  the  uniaxial  cyclic  curve  by  a 
Ramborg-Osgood  relation.  Figure  (3a)  compares  the  test  results 
with  both  the  reference  test  predictions  and  the  exact  numerical 
procedure.  Very  close  agreement  is  obtained.  Figure  (3b)  shows 
similar  calculations  for  the  cyclic  twist  of  o  solid  circular  bar. 

The  third  set  of  data  are  conducted  on  simply  supported 
304  stainless  steel  beams  subjected  to  repeated  reversals  of 
central  load  and  central  deflections  at  1100  °F.  Ten  cycles 
of  load  (or  deflection)  are  completed.  The  experimental  results 
are  shown  in  Fig.  (4a,b) .  The  reference  test  predictions  are 
based  on  the  10th  cycle  cyclic  curve  obtained  by  connecting  the 
cusps  of  three  hysteresis  loops  determined  at  +  E  ■»  0.4,  0.6  and 

1%  .  The  cyclic  curve  is  fitted  to  a  bilinear  relation?  £  s 

+  4  (<r  -  o^)J  /E  where  v  is  a  material  yield  stress  “10000psi 
C  9  1  .  Using  this  relation,  the  reference  strain  predictions  for 
the  central  load  P  after  10  cycles  of  reversed  deflection  S 
■  i  0.1"  in  a  simple  beam  with  b  -  1",  h  *  2"  and  2L  »  24" 
gives  5290  lbs  as  compared  with  5480  lbs  from  the  test  results. 

An  excellent  agreement  of  less  than  4%  error.  In  the  second 
test,  a  simply  supported  beam  is  subjected  to  ten  repeated 
reversals  of  central  load  p  ■  +  2000  lbs.  The  references  stress 
prediction  for  the  10th  cycle  deflection  range  gives  0.082"  as 
compared  with  an  experimental  deflection  of  0.113".  A  less 
satisfactory  agreement  than  the  deflection  controlled  test  but 
still  may  be  accepted  in  view  of  the  smaller  plastic  strains 


involved  in  the  load  controlled  test.  It  is  interesting  to 
note  that  the  present  predictions  compare  very  well  with  a 
full  transient  analysis  of  the  beam  problem  [9]  . 

A.  Cyclic  Creep  Analysis 

The  phenomenon  of  cyclic  creep  is  the  cyclic  strain 
growth  observed  in  metals  when  subjected  to  the  combined 
effect  of  steady  and  cyclic  load  components.  A  closely 
related  phenomenon  to  cyclic  relaxation  of  stress  under  the 
combined  effect  of  steady  and  cyclic  strain  components. 

Plots  of  cyclic  creep  strai  against  cycles  are  quite 
similar  to  the  conventional  time-dependent  creep  curve 
with  a  primary  stage  of  decreasing  strain  rate,  a  secondary 
stage  of  constant  strain  rate  and  a  tertiary  stage  with 
increasing  rate  before  the  specimen  fails  due  to  fatigue. 

The  strain  growth  is  usually  accompanied  by  a  hysteresis 
loop.  Hence,  the  total  plastic  strain  during  a  typical 
load  cycle  is  composed  from  a  small  reversed  component  and 
a  drifting  growth  component  [4]  . 

The  steady  cyclic  rate  of  strain  growth  depends  on  the 

* 

maximum  stress  and  the  stress  amplitude  flir  of  the  stress 
cycle  1 10]  ,  i.e. 

CN  =  A  tr*n  .  Atrn^  (14) 


Jests  carried  out  on  different  grades  of  commercially  pure  Copper  1 10 , 11] 
show  that  the  values  of  n,#  depend  largely  on  the  level  of  tr*  as 
compared  with  the  material  yield  stress.  For  example;  n=3.2, 
X  =  0.03  at  cr*/®y  =  (0.5-1),  n  =  4.6,  *  =  0.2  at  r*/°y  =1-15 
and  ns18,  X  *03  at  A  where  o-y  i  s  the  1%  yield  stress.  Hence, 

it  may  be  generally  stated  that  for  the  range  of  maximum 
stress  encountered  in  design,  the  values  of  24  n4  10  and 
O.024  ^40.2  are  acceptable. 

The  material  law  in  eq.  (14)  is  used  here  to  estimate 
the  effects  of  cyclic  creep  on  deformations  of  structures 
subjected  to  cyclic  loading.  The  solution  is  expressed 
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in .  a  reference  stress  form.  This  is  achieved  throuqh 

using  the  so-called  rapid  cycle  solution  applied  successfully 
in  the  analysis  of  time-dependent  creep  deformation  in  com¬ 
ponents  subjected  to  cyclic  loads  with  short  cycle  times 
compared  with  the  material  time  scale  [12J .  The  technique 
provides  an  upper  bound  on  cyclic  structural  deformation. 

In  rapid  cycle  solution,  the  cyclic  stress  field  in  the 
structure  is  assumed  to  be  composed  from  the  cyclic  linear 
elastic  stress  field  and  a  residual  stress  field  which  is 
in  equilibrium  with  zero  external  forces.  The  magnitude  of 
the  residual  stress  field  is  determined  by  enforcing  the 
compatibility  of  cyclic  strains  resulting  from  the  cyclic 
stress  field  over  a  complete  cycle  of  load.  The  theoretical 
basis  for  the  rapid  cycle  solution  may  be  found  in  ref.  fl2] 
and  its  application  to  the  analysis  of  cyclic  growth  of 
structures  are  found  in  ref.  [10,  11]  . 


As  an  example,  consider  a  rectangular  beam  subjected  to 
the  cyclic  history  of  bending  moment  which  varies  between  M1 
The  cyclic  elastic  stress  field  is: 


<r(y)  =  |  — ~3  .  y 
bh 

..  M_ 
a  #.  3  2 

<r(y)  -  2  — t  •  y 

bhJ 

3  MrM2 

A«r(y)  =  4  -^4  .y 

bhJ 


during  1st  half  cycle 

during  2nd  half  cycle 


(15) 


A* 


Hence  the  cyclic  stress  field  a-  is  assumed  as  the  sum  tr  +  j5  i.e., 


3  M1 

cr  ly)  =  *•  — *-y  +  _p  (y)  during  1st  half  cycle 
bh 

cr  (y)  =  j  - _.y  +.P  (y)  during  2nd  half  cycle 

*  bh 


(16) 


where  P  (y)  is  the  residual  stress  field  which  should  satisfy 
the  equilibrium  conditions 


J"cr*(y)  .  dy  =  0 
_h 


(17) 
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Compatibility  of  cyclic  creep  strains  resulting  from  tr  (y)  and 
o*  (y)  requires  that; 

-,n  r- 


•5 

l  — ^  .y+p(y) 
bh 


3  MrM2 


bh' 


3n 


"  KN^ 


(18) 


where  is  the  cyclic  curvature  rate .  Hence, j3  (y)  is  given 

by  : 


?  (y)  = 


V 

1/n 

' 

1 

A 

3 

M1-M2 

It' 

bh3 

( l-£n) /n  3  M1 
•  y  o  T*  y 

*  bh 


(19) 


Upon  substitution  in  the  equilibrium  conditions,  Kjj  is  obtained 
as ; 


kn  - 


=  A 


—  - 

n 

r  _ 

n 

r  -r 

2n+l-n# 

1 

3M^ 

3(M1-M2) 

3n 

h  * 

2bh2_ 

4bh2 

L_  -J 

(20) 


It  is  noted  that  the  associated  cyclic  stress  field  may 
possess  a  singularity  at  y=0  which  can  be  seen  from  eq.  (19) 
for  n!(Al.  The  singularity  may  be  overcome  by  introducing  an 
assumption  that  the  maximum  stress  at  y=0  does  not  exceed  an 
upper  limit  which  may  be  taken  as  some  multiples  of  the 
material  yield  stress.  A  similar  procedure  has  been  used  in 
ref.  110]  in  the  analysis  of  ratchet  strain  rates  in  a  circular 
Bree  element. 

The  curvature  relation,  eq.  (20),  may  be  expressed  as: 


K.,  =  A(<VerK  )n  («C  Acr  )n* .  i  f(n,$) 
N  1  max  2  max  h 


where 


max 


-  2  Vbh  '  -  4  (M1~M2*/bh  • 


£(n,8) 


2n+l-n2( 


3n  0<l'<<2  J 


(21) 


(22) 


The  values  of  £>h 


can  bo  chosen  for-  the  oract: 


ranges  of  2  Lr.L.  10  0.02  L.  X  j .  2  such  tn?  ..  r 

is  weakly  dependant  upon  n  and  X.  Taking  f  1.2,0. 02)  -  f  ( 1  o  ,  C  .  1 . : ; 
and  f(2,0.2)  =  f  (10, 0.2)  prov'd*'  r-.-c  or; -•>  :o  zf.  and 

c<2  which  yield  their  respective  values  as: 

ot.  -  0  r9  2  ,  oc  =  0.671.  The  corresponding  variation  of  f 

X  4. 

with  n,X  does  not  vary  from  the  mean  value  -•  1.506  by  more 
than  +  2.8%  as  shown  in  table  M)  .  Hence  one  final  result  ioi 
reference  stress  solution  is: 


=  1.596 


K  Cn[°- 


671  cr  ,  0.592  Ao- 

max  max 


(23) 


Table  (4)  Values  of  f(n,#)  used  in  cyclic  creep  of  bear..; 


*  *=  0.02 

/  ._  o  > 

J  0 

n=2 

1.552 

1.58 1 

1.612 

n=4 

1.575 

1. 60S 

,r  /.  r 

n=6 

1.573 

1.611 

1.638 

00 

11 

C 

1.564 

i .  5  0  4 

r.627 

n=10 

1.522 

1, 595 

X  -  u  i  .' 

(  the 

cioseness  of  coefficients 

it;  eg.  (I 

!  3 )  the; 

in  the  case  of  cyclic  plastic*. tv;  Table  (3) 


those  obtained 


An  analogous  result  for  cyclic  creep  deflection  •>_'  n  cantilever 
under  an  end  load  which  caries  between  p,  and  p  can  be  easily 

j  £ 

obtained  as: 


SN  =  0.466  £-•  En 


r  p ,  ■ 

2  b  h‘7  3 


■j .  a 


I  if ' 

7  1  —i- 


J-  u 


?  , 


I  2  ■ 


-j .  Conclusions 

The  reference  stress  method  can  be  apn.lj.ed  to  anal"- i-  or 
structures made  from  cyclically  hardened  materials.  The  technique 
is  extremely  useful  during  the  initial  staor-n  of  design  sir.ee 
reduces  the  problem  to  conductin',  a  .  *  il^ 
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hardening  test  at  the  reference  value  of  stress  (or  strain) 
amplitude.  The  accuracy  of  the  technique  has  been  demons¬ 
trated  through  comparisons  with  published  test  results  and 
exact  numerical  solutions  for  simple  beams.  Checks  on  more 
complex  components  such  as  plates,  panels  and  tubes  should 
be  conducted. 

The  applicability  of  the  RSM  to  analysis  of  cyclic  creep 
behaviour  of  structures  is  demonstrated  by  extending  the  use  of 
the  so  called  rapid  cycle  solution  used  in  time-dependent  creep 
analysis.  Tests  on  components  operating  under  conditions  where 
cyclic  creep  is  dominant  are  required.  In  principle,  the  RSM 
can  be  applied  to  analysis  of  structures  under  cyclic  relaxa¬ 
tion  conditions. 
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